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1. Introduction 


In the conventional approach to non equilibrium thermodynamics, as _pre- 
sented for instance in Casimir’s review article [Casimir, 1945] or in de Groot’s 
monography [de Groot, 1951], based on Onsager’s classical work [Onsager, 1931], 
the following assumptions are believed to be fundamental: 


. Time reversal symmetry (principle of microscopic reversibility) 

. Fluctuation theory 

The assumption that the phenomenological macroscopic laws of motion (as 
for example the relation between temperature gradient and flow of heat or 
between concentration gradient and flow of substance) hold also in the range 
of fluctuations. 


a 


We shall examine these conceptions very briefly. 

The principle of microscopic reversibility expresses the invariance of the fun- 
damental microscopic equations of motion of a closed system under time rever- 
sal. It holds in quantum mechanics as well as in classical mechanics. The con- 
cept of time reversal symmetry was probably first used in statistical mechanics 
by Loschmidt, who used it in a negative sense as an argument against Boltz- 
mann’s H-theorem. Einstein used detailed balance arguments, whose origin is 
time reversal symmetry, in various connections. Already the early discussions 
showed that the application of time reversal symmetry to irreversible processes 
is not completely trivial. A fundamental difficulty is the interpretation of mac- 
roscopic measurements on a microscopic level. This problem can be formally— 
and is usually—neglected as far as classical mechanics can be used for the mi- 
croscopic description of the system.! The reason for this is of course that classi- 
cal mechanics does not forbid us to assume that an observation can be made 
arbitrarily detailed and that nevertheless the perturbation of the system caused 
by the measurement can be made arbitrarily small. This means that the system 
can be considered as closed in spite of the fact that observations are performed 
on it. But according to the quantum theory of measurement a system on which 
observations are made cannot even approximately be considered as closed. On 
the other hand time reversal symmetry is strictly valid only if the system is 


1 Interesting exceptions are the papers [Green, 1952] and [Yamamoto, 1953]. A forrunner of 
these is [Zermelo, 1900]. 
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closed. Obviously there is some approximation involved in the application of 
time reversal symmetry to irreversible processes. 

Fluctuation theory is in general used in the form given by Smoluchowski and 
Einstein [Smoluchowski, 1908], [Einstein, 1910], in which the Boltzmann relation 
between entropy and probability is reversed and the empirically known dependence 
of the entropy on observable quantities is used to determine the probability of 
a deviation from equilibrium in an adiabatically insulated system. In the actual 
applications Boltzmann’s relation is used only to compute the statistical second 
order moments, the so called covariances of the fluctuating quantities. 

It should be remembered that the object of Onsager’s original work (in con- 
trast to the more recent papers [Onsager and Machlup, 1953]) was not a theory 
of the kinetics of fluctuations but the proof of reciprocal relations for the coef- 
ficients occuring in the macroscopic equations of flow. The concrete problems 
considered by Onsager were in fact: Heat conduction in anisotropic crystals, and 
chemical monomolecular triangle reactions. Fluctuation theory was used only as 
an, ingenious, artifice for the purpose of obtaining reciprocal relations. 

It was particularly stressed by Casimir that Onsager’s assumption that mac- 
roscopic equations applying originally to deviations large compared with fluctua- 
tions, can also be applied to the average behaviour of these fluctuations them- 
selves, is really a new hypothesis. Casimir tried to justify the Onsager assumption 
qualitatively by associating with the fluctuations a characteristic time interval 
which is required to establish a state of steady flow. This time interval is as- 
sumed to be small compared with the time in which the deviations from equi- 
librium is appreciably reduced. That means that the macroscopic equations are 
assumed to be valid for the average behaviour of fluctuations on a ‘‘coarse 
grained” time scale. What Casimir had in mind was apparently an approach to 
irreversible processes from the kinetic theory of gases by use of the Chapman- 
Enskog method. Perhaps the most instructive illustration of the necessity of dif- 
ferent time scales is however contained already in the beautiful work of Smolu- 
chowski on the connection between diffusion, Brownian motion and coagulation in 
colloids, [Smoluchowski, 1916]. 

At present the attempts to establish a direct contact between the kinetics of 
fluctuations in equilibrium systems and the macroscopic equations of flow in non 
equilibrium systems are of a rather tentative nature. It therefore seems sound 
to attack these two problems independently. 

Whereas a great number of papers have been concerned with the develop- 
ment of fluctuation theory on the basis of Onsager’s hypothesis,! the possibility 
of a direct approach to the macroscopic equations of flow has been almost 
completely neglected, if we leave out of consideration those works which ap- 
proach the problem through kinetic theory (Boltzmann’s equation) or make use 
of a master equation and thus replace the original problem by a new, as yet 
unsolved one, the derivation of the master equation. An exception is the work 
of Bergmann and Lebowitz [Bergmann and Lebowitz, 1955], [Lebowitz and Berg- 
mann, 1957]. These authors consider a system interacting with a set of infinite 
heat reservoirs—in the second paper also more general reservoirs are considered. 
The internal dynamics of the system is assumed to be described completely by 
its Hamilton function. Its interaction with the reservoirs consists of random im- 


' In particular [Onsager and Machlup, 1953] and [Hashitsume, 1952, 1956]. 
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pulsive ‘“‘collisions’’. Mathematically this is stated in terms of a generalized Liou- 
ville equation for the propability density in phase space containing an integral 
operator (collision term) which takes account of the random interactions with 
the reservoirs. 

The first generalization of Liouville’s equation to include the action of a heat 
reservoir on a system was given by Klein [Klein, 1921, 1922]. In his theory the 
effect of the random interaction with the surrounding is taken care of by the 
introduction of a generalized friction and diffusion term into Liouville’s equation. 
Klein’s equation may thus also be considered as a generalization of the Smo- 
luchowski-Einstein equation for Brownian motion. Klein’s equation generates the 
transition probability density of a Markov process. The corresponding stationary 
Markov process is symmetric with respect to time reversal, and may be used 
as the basis of a kinetic theory of fluctuations. : 

In Klein’s as well as in Bergmann’s and Lebowitz’ theory the interaction term 
is constructed so that the probability distribution of the system will in the pres- 
ence of one reservoir approach a Gibb’s canonical distribution and in the pres- 
ence of several reservoirs approach a stationary distribution, which is in general 
not canonical. It can then be demonstrated that Onsager’s reciprocal relations 
are satisfied by the coefficients for heat conduction if the temperature differences 
between the reservoirs are small. We shall show in this paper that the particular 
model used by Bergmann and Lebowitz is not necessary for the derivation of 
Onsager’s relations. In fact no special model is needed at all. This is of some 
importance as the Bergmann Lebowitz formalism is not easily extended to quan- 
tum mechanics! and also because objections can be raised against the introduc- 
tion of a phenomenological term into the Liouville equation as this implies a mix- 
ing of two different time scales, one fine grained corresponding to the internal 
dynamics described by the Hamilton function, and one coarse grained corresponding 
to the interaction with the reservoir.” 


2. Outline of some problems connected with Onsager’s reciprocal relations 


In this section we shall give a survey of some typical fluctuation and tran- 
sport problems which will be closer examined in later sections from the point 
of view of general statistical mechanics. 

We consider a thermodynamic system which is originally completely closed. 
Let the macroscopic state of the system be characterized by the values of a set 


of extensive quantities A’, ..., A” 
Sen gp A a (1) 


Illustrative examples of the extensive quantities A are internal energy and vol- 
ume of subsystems, number of moles of different species in subsystems etc. The 
system being closed, the values of the A:s are fixed and show no fluctuation 


(by definition). 


Fluctuations in an open system 


Now let the system be brought into contact with a reservoir with which it can 
exchange the extensive quantities corresponding to A. The reservoir is assumed 


1 See [Willis, 1957]. .— 
2 Of course such a mixing occurs also in the Boltzmann equation of the kinetic theory of gases. 
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to be so large compared with the system that its macroscopic state is not affected 
by this interaction. After some time the system reaches a macroscopic equilibrium 
state which is characterized by the values of the intensive quantities Baer 


B, = Bi, tee Br =Bm (2) 


which correspond to the extensive quantities 4 and determine the macroscopic 
state of the reservoir. 

Let the values of the quantities A be measured at different times, after that 
equilibrium has been reached, and denote by 


(G5) weet) =a) (3) 


the values observed at the time ¢. In general the values a (t,) and a (t,) obtained 
at different times ¢,<¢,, are different and it is natural to call the difference 
a(t,)—a(t,) the flow of A-quantity from the reservoir into the system during 
the time interval (¢,, ,). However when taking a very large number of obser- 
vations we will find a certain stationary behaviour of the observed values. This 
is precisely the criterion for equilibrium. Assuming that the values a (t) are pro- 
duced by some random mechanism we may try to construct models for this 
mechanism. Ideally such a model should be derived from the microscopic equa- 
tions of motion, but except for very special cases this has not been possible 
so far. On the other hand it is usually easy to express invariance properties 
of the microscopic equations as symmetries of the random process associated with 
the observation series. In particular this is true for time reversal invariance. 


Reaction of the system to a sudden change of the constraints 


The contact between the system considered above and the reservoir is broken 
and the system is brought into contact with a second reservoir of the same 
kind as the first one but in general with different values of the intensive parame- 
ters B. Originally the system is not in equilibrium with the second reservoir. As a 
consequence it will be possible to registrate average flows of the quantities A 
from the reservoir into the system. It is quite natural to consider the differences 
in the intensive parameters B as the “driving forces’ producing these average 
flows. The dependence of the average flows on the differences of intensive para- 
meters is restricted by the time reversal symmetry (valid for the closed system) 
and for small differences, that is in the domain of a linear approximation On- 
sager’s relations are valid. This will be demonstrated in this paper. 


Steady state in contact with several reservoirs 


The system considered above is brought into contact with several reservoirs, 
all of the same kind, but with different values of the intensive parameters B. 
A steady state, which is in general not an eqilibrium state, is assumed to be pro- 
duced in the system in the course of time. In this steady state the average 
net flows of the quantities A from all the reservoirs into the system are zero. 
This means that the different reservoirs are connected by steady flows. The dif- 
ferences of the intensive parameters B may be considered as the “driving for- 
ces” of these steady flows. The dependence of the average flows on the differ- 
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ences is restricted by the time reversal symmetry and in the domain of a linear 
approximation Onsager’s relations are valid. The demonstration is given in 
this paper. 


3. Statistical description of an open system 


We shall consider the experimental situation described above as fluctuations 
in an open system. The intensive parameters B which characterize the state of 
the reservoir are chosen as the conjugated variables of the extensive variables 
A with the entropy function, or possibly a generalized entropy function, used 
as thermodynamic potential. If for instance the A-variables are 

2 ore 71) 1) () 

PES DO ee aN eae wera at eel etl ot ats vs (4) 
that is internal energy of subsystem (1), ..., volume of subsystem (1), ... , number 
of moles of the substance 1 in subsystem (1) ..., the corresponding B-varia- 
bles are 


Here 7’ denotes temperature, p pressure, and yj chemical potential. Upper in- 
dices refer to subsystems and lower indices to species of substance. As indicated 
by the notation we do not necessarily assume that different subsystems are kept 
at the same temperature, pressure, or chemical potential. 

We shall assume that the state of the system at any time is determined by 
the measurement of an observable Z={Z(1), ..., Z(n)} whose components Z (j) 
take values on the real axis. The extensive variables A are assumed to be unique 
functions of the observable Z 


ie) 5a O (Coen (Cin GNC Gn (6) 


The observed value of Z(j) as a function of the time may be considered as 
a realization of a stochastic process {z(t,j);— co <t< +0, 7=1,2, ..., n}. It is 
mathematically convenient and physically irrelevant to have the stochastic process 
defined for all values of the time parameter ¢ on the real axis. 

The stochastic process z(.,.) associated with the observation of Z is deter- 
mined by the family of all finite dimensional distribution functions of the form 


0 (ts Iu Ste a Ps ate LE) (7) 


representing the probability of finding z(t,7,) in the interval J, of the real 
axis, ..., and z(t,,4,) in the interval J;,. 

If the system has been in contact with the reservoir for a very long time 
before the measurement of Z starts we are entitled to consider the stochastic 
process z(.,.) corresponding to the observed values of Z as stationary. That is, 


we assume that 
Olt, at ba Iys eee yg babe Nes 1 coe L,)=0 (tp eee yg bres Iie3 Le eve I,) (8) 
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for every finite dimensional distribution function and every real t. In particular 
it follows from (8) that the distribution function corresponding to a single instant 
is independent of the time 


ONS eae ire oul Penn 1 onl ere PEG (9) 


and that the distribution function corresponding to two instants depends on 
their difference only 


OG Wyle nla tay Lr ng ths lyn ean) Oa lg te ee 5 lon) (10) 


We shall assume that each component of Z has a definite parity under time- 
reversal (reversal of the direction of the microscopic motion): 


0 ()= + 1 or 1. (11) 


If the composite system, including the reservoir, is a closed dynamical system 
described by classical mechanics or by quantum mechanics and if the stochastic 
process z(.,.) is stationary one can prove that the stochastic process will also 
possess the symmetry 


C= tts eee 5 Steins Oi) tae eee 6 (ji) Le) = 0 (ts hh RIEFOR thes Jes nee eee yg I,,). (12) 


The result of 6(j) applied to an interval J=(a,b) is J if 0(j)= +1 and (—6, 
—a) if 6(j7)= —1. For the particular distributions (9) and (10 we obtain in the 
case of time reversal symmetry 


DO, (ON) da con gO (A) iy laaeeee gas (13) 


showing that the stationary (equilibrium) distribution is invariant under time 
reversal, and 


Gg =EO UD) Lats Or) LOL) [nai oy OL eno, G01 egeenan 
O(a) Leg si CL) Le src veg ON 1) Legs) Sag bod ae) eacichl el made eee ela 


The proof of the symmetry relation (12) will not be given in this paper. 

The values /, ..., 8, of the intensive variables B will appear as parameters 
in every finite dimensional distribution function (7). In particular we shall assume 
that the stationary distribution 9,(.) is of a generalized canonical form [Uhl- 
horn, 1960 a (7)] 


a Leas : 
0, (dd, ... d¢,)=exp fe - Bua (2) CREP Cre ven (15) 


where H is Boltzmann’s constant, and Y=Y(f) is the Massieu-Planck function 
determined by the normalization of (15) and related to the entropy function S («) 
by the Legendre transformation 


Siajse"8,- VB), ea. (16) 
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The properties of the distribution (15) have been discussed at length in [Uhl- 
horn, 1960 a], and we shall make the same assumptions about the functions 
a"(f) as in that paper. The measure yu (d£) occuring in the definition (15) is 
independent of the parameters fj, ..., Bm. In fact w(d¢) is determined either 
by the linear operator in state vector space or by the phase function which is 
associated with the observable Z, that is (dé) is the structure function of Z. 

For the purpose of deriving Onsager’s reciprocal relations we do not need the 
explicit form of any of the distributions (7) except 9, (.). The only properties 
of the distributions which are needed are the symmetry relations (10) and (U2). 


4. Hilbert space treatment of the stochastic process 


In this section the analytical methods needed for the derivation of Onsager’s 
relation are developed. 

Let t=(t,7) denote the parameter of the stochastic process z(.,.) associated 
with the observable Z. We note that o(t, ..., t%;d¢, ... df,) is a probability 
measure in the k-dimensional Euclidean space R*. The Hilbert space of all com- 
plex functions of k real variables whose square is integrable over R* with the 
measure o(t1, -..; T,3 @C) will be called L(t, ...., 1). 

By H(t, ..., tT) we shall denote the Hilbert space of all random variables 
of the form 


X (2 (T), Oe} 2(Tx)); X (¢4, eee 9 Cy) €.D? (T1, Oa) Te) (17) 


with the scalar product 


CK (t,), «0 2G) 1X5 2G), --, 2 = 
= Ae Pee Gg aoe PT is sg Ure Ca keeca ou) 2 (Ls) 
Rk 


and the norm 
|| X (z(t), eae z(t;))]| > CX (2 (71), ree @ (t,,)) |X (z (7), Cut} z(ty))>*. (19) 


It is important to observe that though the Hilbert spaces TG we Tk) and 
H(t, ...,; T;) are isometric, according to the definition (18), they are not iden- 
tical. 

The union of all Hilbert spaces H(t, ....t,) determined by a finite number 
of parameter values of the stochastic process is clearly a linear space with a 
scalar product corresponding to (18), and after completion a Hilbert space Jal fie 
obtained. We note that though the elements of a Hilbert space of the type 
H(t, ..., T) can be analytically characterized with the aid of functions of a 
finite number of real variables this is no longer the case with the elements of H. 
The elements of H will nevertheless be called random variables.' 

We can construct a Hilbert subspace of H corresponding to any set Jot 
parameter values of the stochastic process z(., .) by forming the union of all 
Hilbert spaces determined by a finite number of parameter values which belong 
to the set J. After completion we obtain the Hilbert space H (J) H. 


1 Tt is not difficult to show that the elements of H are random variables in the usual sense of 
the word, that is that they can be represented as measurable functions on a probability space. 
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If « is an arbitrary element in H and J is a set of parameter values of the 
stochastic process z (.,.) the orthogonal projection of « into H (J) is called the 
conditional exceptation of the random variable 2, given {z(t); 7 ET}. For, denoting 
by E’ the projection on H (J) in H, it is evident that E’ x is the best approxi- 
mation of « by an element in H(J) if the goodness of the approximation is 
determined by variance of the random variable x— 7x. That is Hx is the least 
square estimate of 2. In particular the unconditional expectation of x is equal 
to its projection onto the Hilbert subspace of constant random variables. This 
subspace is one-dimensional and is spanned by the unit vector e=1 which 
evidently is contained in all Hilbert spaces determined by a finite number of 
parameter values. Hence the expectation of the random variable x can be 
written 


Hu=ece 


\ 


By (20) 


Sets of parameter values of the stochastic process z(.,.) which are deter- 
mined by conditions on the time variable t only are of particular interest to us. 
By H(t) we denote the Hilbert subspace of H corresponding to the values 
(t,1) ..., (t,m) of the parameter and by L*(t) the corresponding function space. 
Analogously we introduce H (t,,..., t,) and LD? (t,, ..., t,). 

The product of the projections E%'"™ and Ht 


BO ® EG: tH) (21) 


by the isometry (18) determines a transformation from the function space 
LE? (ti, ..., Ti) into the function space L?(t,, ..., 7) which we shall denote by 


Clie ne) Gee (22) 


Being isometrically equivalent to the product of two projections the operator 
(22) is evidently bounded by 1 and in fact its norm is equal to 1 as the unit 
vector e (20) is transformed into itself by the application of (21). Further the 
operator (22) has the easily verified property of transforming a non negative 
function into a non negative function. 

So far the discussion applies to any type of stochastic process. We shall now 
take account of the symmetries (8) and (12), stationarity and invariance under 
time reversal. 

It follows from (8) and the definition of LZ? (t,,7,, ..., t,,4,) that 


TAGS If Site EG =i (CR Aiea (23) 

Hence it is possible to define a transformation U (t) by 
U (t) X (2 (ty91), «++ > & (tes Je) =X (Z(t, +6,9,), ---, 2 (Ge +t, 4,)) (24) 
where X (¢,, ..., ¢,) is any function belonging to the space (23). The transforma- 


tion U (t), being bounded in its domain of definition, can be extended to the 


whole Hilbert space H and gives rise to a one-parameter group of unitary 
transformations 
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U (yr 8) (jr = Oi) 80) (t,) U (ty) = U (t, + ty). (25) 

We shall make the additional assumption that if 2 is an arbitrary element in H 
the expression <x|U (t)|x> is a continuous function of f. 

In order to take account of the time reversal symmetry (12) we define the 


operator 7’ by 


PX (Gr 95 anaes ben) ee aH) On et), ree ON Gs.) = tau) (26) 


for any function X (¢,, ..., €,) belonging to the space (23). The definition is 
consistent as it follows from (12) that the function X (6 (j,) C4, ..., 6(j.)) be- 
longs to L?(—#,,4), ..., —t,,j,). Being evidently bounded on its domain of de- 


finition 7’ can be extended to H. It follows from (12) that 7’ is unitary 
(ap t=, (27) 


The involutory property of 7 is a direct consequence of the definition (26). 
From the definitions (24) and (26) we obtain the connection between the 
group of “kinetic transformations” U (ft) and the time reversal operator T 


TU (t)=U(—2)T. (28) 


We shall conclude this section by a study of the operator C'(t,,t.) which is a 
particular case of the operator (22), defined by (21). The stochastic process z(., .) 
being stationary, the function spaces L?(#,) and L?(t,) are identical both corre- 
sponding to the distribution function 9,(d¢) defined by (8), and denoting this 
space by Li, we find that C(t,,t,) is a linear transformation of L{ into itself. 
From the unitarity of U(t) it follows that 


U (t) B® U (t)* = BO, (29) 


and hence the operator C(t,; t,) can depend only on the difference t, —4,, so that 
we may write 


C(t; te) = C(t, — ty). (30) 
Recalling the definition of C(t,; t,) we find that C(t) has the properties 
C(0)=1; C@)*=C(—2). (31) 


The invariance of the stationary distribution g,(d¢) under time reversal (13) 
makes possible the definition of the operator 7’ in ik 


TK (Cy ++, En) =X (O(a, «-- » O(m) Cn). (32) 


The time reversal operator 7 in L{ is obviously involutory and unitary. It is 
connected with C(t) by the relation 


POKL = C(t), (33) 
which is easily obtained by applying 7 (in H) to the projection B® and consid- 
ering the definition of C (é). 
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The physical meaning of the operator C'(t) in L* is the following. Let the value 
of the observable Z be observed at the time f) and let X (Z) be an arbitrary 
complex function of the observable Z belonging to the space ie We denote the 
observed value of Z by ¢. The expected value of the observable X (Z) at a time 
t) +t is then given by 


Cw) X (0). (34) 
We note that with X(f) and Y(¢) in Li we have identically 


<X (2 (t1)) | Y (@ ())> =<X @ (4) |U (2-4) | YG) = 


=(X()O(t-t4) Y(0)e, (dL). (35) 


Rn 
As we have assumed U(t) to be (weakly) continuous the same holds for C(t). 
In particular 
lim C(¢)=1. (36) 


t 0 


There is no restriction from the physical point of view in assuming the sto- 
chastic process determined by (7) to be ergodic, or even mixing (compare the 
discussion given in [Uhlhorn, 1960b]). Hence the behaviour of C(t) for large values 
of ¢ is given by 


lim C(t) X=[X(C)o,(46)=EX, (37) 


t>+0 


expressing the approach to equilibrium (decay of fluctuations). 


5. Fluxes 


It was assumed in section 3. that the extensive variables A are unique func- 
tions of the observable Z, equation (6). Having associated the observable Z with 
a stochastic process corresponding to observations of Z when the system is in 
contact with a specified reservoir we obtain through the functional relation (6) 
a stochastic process associated with the extensive variables 4 


G2 (6, AG (2b yp me gO BE, 0) hs (38) 


The derived stochastic process (38) is obviously stationary with z(.,.) and pos- 
sesses the time reversal symmetry. In fact the last symmetry relation is satisfied 
by the process (38) in the strengthened form of detailed balance. For it follows 
from the form of the equilibrium distribution 9, (dé) given by equation (15) that 


the functions a () must be invariant under the time reversal transformation (32) 
The expression 


[al (2(t-+8, J) — a" (2, .))](s > 0) (39) 


is the net flow per unit time of A"-quantity from the reservoir into the system 
during the time interval (¢,t+.s). The mean value of (39) is zero as a consequence 
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of the stationarity of the stochastic process. But the conditional mean value of 
(39) given that 

z(t, .)=C (40) 


will in general be different from zero. Using the conditioning operator C(t) de- 
fined by (30) we can write the conditional mean value of (39) as 


1 
5 Le(s)— Na*(¢). (41) 


This does not necessarily require that the function a“(.) belongs to Li as it may 
be possible to extend the domain of definition of C(t) beyond this class. 

The flux of A“-quantity from the reservoir into the system will be defined as 
the stochastic process J“(z(f,.)), where the function J“(¢) is defined by 


A aeatnes es Nae (Cas (42) 
syo § 


The limit (42) is assumed to exist in the sense that for arbitrary X(.) in L} 


Ieee Fe z ‘ 
lim —<X|C(s)—1|e =<X| J*. (43) 
syo § 
It should be noted that J”(.) depends on the stochastic process z(.,.) through 
C(.), and hence depends on the particular reservoir. As already mentioned the 
expectation of J“(z(t,.)) is zero 


EJ" (z(t, .))=0. (44) 


Applying the relation (33), which expresses the time reversal symmetry of the 
operator C(.), to the definition (43) of J”(.) and taking into account that 


Lo =a. (45) 


we find 


j= = ie, (46) 


that is the fluxes are odd with respect to time reversal as is to be expected. 

The discussion given above applies when the system is in contact with any 
reservoir of the same kind. We shall distinguish between quantities referring to 
different reservoirs by an index. For instance J() (.) is the flux defined by the 
stochastic process z(.,.)qa) corresponding to contact with reservoir number (1), 
with the parameter values f1,1), ... 1 Dams: 


6. Reaction of the system to a sudden change of the constraints 


The physical situation has been described in section 2. The system has been 
for a long time in contact with reservoir number (1). At a certain time ton 
is replaced by another reservoir of the same kind (2). The observations of the 
“state” of the system consist in the measurement of the observable Z (section 
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3.), which is associated with the stochastic processes 2(., .)a) and 2(.,.)2) when 
the system is in contact with reservoirs (1) and (2) respectively. For times prior 
to f) we may use the stochastic process 2(., .)a) to describe the observed values 
of Z. The stationary distribution 01,)(.) is the generalized canonical distribution 
(15) determined by the parameter values fia), ... Ama. After the instant t) the 
observed values of Z cannot in general be associated with a stationary sto- 
chastic process. It is however reasonable to assume that it is possible to describe 
the situation after t,; by the conditional process derived from the stationary sto- 
chastic process 2(.,.)2). As the last process is assumed to be mixing in the sense 
of [Hopf, 1937], which in particular implies (37), the system will ultimately reach 
a new equilibrium state corresponding to a generalized canonical equilibrium 
distribution 01)(.) determined by the parameter values fi, ... ,Bm@). 

As a consequence of the non stationary character of the stochastic process 
associated with the observation of Z, the mean values «', ...,«” of the stochastic 
variables a'(t), ...,a”(t), defined by (38), will be time dependent for t>t). Only 
if the two sets of parameters fa) and fiz) are equal will the mean values be time 
independent. We can therefore expect a phenomenological relation of the form 


at! (to +t) — a (to) = LIP (t) AB, (A Bu = Bue Bua)? (47) 


to be valid to the first order in Af,. 
We shall show that the coefficients 


Ly ={ : Lato +) —a"(o]| ee) 


obey Onsager’s reciprocal relations 
L” (t) = L(t). (49) 
Proof of reciprocal relations 


Denote the index of probability corresponding to the equilibrium distributions 
ia) and Qi) by 


Gory (C Y (Bo) =x (na (¢)] (n= 102). (50) 


According to our assumption, which is justified by the largeness of the reservoir 
(2) and the smallness of Af,, we have 


a (to +t) — a" (to) = [[C (t)) — la” (2) exp fou (C)} w (d6), (51) 


where C'(t)2) is the conditioning operator of the stochastic process 2(., .)). The 
relation (48) is equivalent to 


py __ { Y lM ] lu 
L l oe Ke aek (aces a) tea (in)]! : (52) 


Bay=By 


1 Summation convention, w=1,...,™. 
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and inserting (51) into this expression we obtain 


ad 


v 1 Y Lys V/s f 
LE ~ | (6 Qe — la" (C)a" (2) exp {oe (£)} u (dC), (53) 


“ 


» 
n 


Bay = Pe implies oy () = 62) (C). (54) 


The expression (53) can also be written as a scalar product in Li) 
Pee 
== <a"|C(@—1|a">e, (55) 


and using the properties (31) and (33) of O(t)s), recalling (45), we find that (55) 
is equal to 


1 Af GX v ve 
is |C (Qe — 1] a> = L. (56) 
This is the desired result (49). 


7. Steady state in contact with several reservoirs 


The physical situation was described in section 2. The system is in simul- 
taneous contact with a number of mutually independent reservoirs (r=1, 2, ...) 
of the same kind but with slightly different values of the parameters B 


Bie, Pao (poo yer a): (57) 


It is assumed that macroscopically the system reaches a steady state. This steady 
state is an equilibrium state only if the parameter sets (57) corresponding to 
different r are equal. 

As above the observation of the system consists in the measurement of the 
observable Z. The series of observations of Z is associated with a stationary 
stochastic process, stationarity of the stochastic process corresponding to the fact 
that macroscopically the system is in a steady state. As however the total sys- 
tem, including the reservoirs, is not in equilibrium we cannot assume that the 
stationary stochastic process possesses time reversal symmetry. We denote the 
stationary stochastic process by z(.,.)9, and write 


Bug = Ppt A Bur (i= ae 00; Pl ae oS) 


We are interested in the case when the differences Af, 7) are small. Let the 
equilibrium distribution of the stochastic process z(.,.)9 be determined by the 
index of probability o(.) which for vanishing A {,,~) reduces to the canonical form 


oo(t)=+ [¥ (8) — Baa (C)]. (59) 


x 
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Further we introduce the stochastic processes z(.,.)@)("=1, 2, ...) which corre- 
spond to the observation of Z when the system is in contact with a single re- 
servoir (r). These processes are stationary and possess time reversal symmetry. 
Their equilibrium distributions are of the generalized canonical form. We introduce 
the derivative of the conditioning operator 


Mey = Ibnan S [C (b)ry aa fe (60) 
tyo ¢ 


This limit is assumed to exist in the sense 


ee oe S 
(X | Mo)| Y= ne i C(Hey-1 | Yay (61) 


for certain X(.) and Y(.) in Liq). It follows at once from the properties (31) 
and (33) of C(t)q) that 
CK Mes | Yoel Y | Mnlt Oe (62) 


if X(.) and Y(.) are real functions. 


The superposition assumption 
We shall assume the following connection between the stochastic processes 


z(.,-)o and 2(., Je) (r=1, 2, ...). In the limit of small differences A Bi), for X(.) 
in the domain of definition of WZ ) 
S [Mon X (6) exp {o(6)} u(a2)=0(| ABI). (63) 


— 


The condition (63) is to be understood in the sense that 


Le MX (£) exp {0 (C)} w(d6)]ap-0=9, (64) 
and a » [ tox© exp {a (C)} w(dC) =0 (65) 
OA Byer r=1 ‘ as ~ JAB=0 


The first relation (64) is a direct consequence of the fact that for A fb=0 the 
index of probability o(.) reduces to the canonical form (59), and the fact that, 
the stochastic processes z(.,.)r, being stationary, the mean value of X (z(t, Jor) 
is independent of f: 


oO . EGE 
BEX lb dn) = lim | | [C Wn X (6) — X(Q)] exp {ow (2)} w (AO) = 


=| Me) X (C0) exp {om (O}u(de)=0. (66) 


The second relation (65) on the other hand represents a new assumption, though 
a rather weak and natural one. Its precise meaning will become clear below. 
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The flux of A”-quantity from the reservoir (7) into the system has been de- 
fined in section 5. as the stochastic variable 


Jor) (z(t, Nos Jr) (¢) = Miya" (Gy: (67) 

The mean value of this stochastic variable in the steady state will be denoted 

Jy) = E Fl (2 (b Jo) = | Mena" (£) exp {o(6)} (42). (68) 

As a consequence of the superposition assumption (63) we obtain the necessary 
condition for a steady state 

SF =O ABP, (69) 

expressing the conservation of A“-quantity. It should be observed that it is only 

in the limit of small differences A£,,;) possible to define the fluxes from indi- 

vidual reservoirs (67). Hence it is sufficient that (69) holds to the second order 


in A Bur). 


In the limit of small Af.) we expect a linear relation of the form 


JG, => Less Apes OAS 


a (70) 
between fluxes and affinities. We shall prove that the phenomenological coeffi- 
cients 
o 
in) = |za I 71 
CFIC ; Bacesy | Ag=o ( ) 
satisfy Onsagers’s reciprocal relations 
Loy 9) = Ley ft): (72) 


Proof of reciprocal relations 
Take X (¢)=1 in (62). In view of (66) 


f Mey ¥ (0) exp {om (6)} w(do) =0. (73) 
Differentiate (73) with respect to By: 


| be Ma r| exp {oir (f)} (df) — 
0 Bur) 


‘* | [Mir ¥ (C)] a* (C) exp {or) (f)} (OT) 


and use the symmetry relation (62), recalling that o@) (€) and hence a’ (C) is 
T-invariant : 
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| be Mav ; exp {o@ (C)} u(aC) = 
= = [Mr Y (C)]a* (C) exp {oa ( C)} fa (dg)= 
NA 
_! Ro VE ae 
= JM a* (C)] exp ou (S)5 4 ( ¢) 
xt a 


hee IG (6) exp {on (2)} uw (a2). (75) 


Take in particular Y (¢)=a’"(¢): 


Caen . js 
| etm exp {oe (C)} w(dl) = 


1 re 7 
=2{ totter exp {or (C)} uy (d¢) 
1 lL e v/s e\) so 
sy Jim (C) a" (C) exp {oir (C)5 uu (df) = 
-| |e Int) oe! (76) 
- 8 By ’ 
Let (0) d o()| (77) 
ev Ounn = 
(Y)I & A) A Bucy & Ap=0 


and consider the superposition assumption (65) 


‘ o 
O= Ss Isla, 1S exp {ia ( u( 2] = 
baa ip fat | oF (6) : : sald Ap=0 
OMw Ye : : 
= i sae font} ed) 4 
0 Bury Ap=0 


‘uz B | Or Fj] Ou (C) exp {our (C )} fad (dc 3] 


. (78) 
Ap=0 
Combine (78) with (75): 
S [I TY (£2) Jt (6) exp {oe (2)} w(dL)]as-0= 
ene [> J Mor Y (6) ou () exp {oe (2)} (do) ]aa-o- (79) 


Take in particular Y(¢)=Z'oq», and use the symmetry relation (62) 


= 


358 


ARKIV FOR FYSIK. Bd 17 nr 19 
[J or (C) Je, (C) exp {om (2)} u(d2)]ap-0 = 
[1 LA! (Moy T oy (6)) Gru (C) exp {our (C)} ue (42)]as-0= 


r 


| (Merry T Oru (¢)) O(r')y (¢) exp {oe (C)} lu (d C)lap=0 = 


a 


a [J Oru (L) Ter (2) exp {oe (2)} w(dL)]ap-o. (80) 


The coefficient L6(,.) is according to (71) 


geese } 
Lye = On | Gos) exp {G(r (C)} - i) at 
O Byer) AB=0 


+ [J Jo) (¢) Ory (¢) EXP {Ou (¢)} lM (dC) ap-o- (81) 


The terms to the right are identical with the expressions (76) and (80) which 
have the desired symmetry, so that finally 


RP Sra ¢ 
Linen = Leryn - (82) 


8. Discussion 


It should be observed that the results of this work have been obtained, 
without any other assumption than the validity of time reversal symmetry for 
observation series in equilibrium systems. This is a reasonable assumption, as 
it can be proved for closed systems in macroscopic equilibrium, and is independ- 
ent of whether quantum or classical mechanics is used for the microscopic de- 
scription of the system. It is evident that any particular model for the interac- 
tion between system and reservoirs, which is in agreement with time reversal 
invariance, will give the results obtained here. It is easy to see that the models 
of Klein and of Bergmann and Lebowitz have this property. 


ACKNOWLEDGEMENT 


I wish to thank Professor Oskar Klein, head of this institute, for his encouragement and 


unfailing interest in the work. 


Institute for theoretical physics, University of Stockholm, Stockholm. January 1960. 


REFERENCES 


BERGMANN, P. G., and Lesowrrz, J. L., 1955, New Approach to Nonequilibrium Processes. 
Phys. Rev. 99, 578. ; oa 

Casturr, H. B.G., 1945, On Onsager’s Principle of Microscopic Reversibility, Rev. Mod. Phys. 
17, 343. al 

Erystern, A., 1910, Theorie der Opalescenz von homogenen Flissigkeiten und Fliissigkeits- 
gemischen in der Nahe des kritischen Zustandes, Annalen der Physik 33, 1275. 


359 


U. UHLHORN, Non equilibrium thermodynamics 


Green, M.S., 1952, Markoff Random Processes and the Statistical Mechanics of Time-depend- 


Company. 

ee cae 1952, 1956, A Statistical Theory of Linear Dissipative Systems. Progr. Th. 
Phys. 8, 461, 15, 369. 

Horr, E., 1937, Ergodentheorie, Springer. 

Kuen, O., 1921, Zur statistichen Theorie der Suspensionen und Lésungen, Arkiv for Matema- 
tik, Astronomi och Fysik 16 N:o 5. 

—— 1922, Zur statistischen Theorie des Temperaturgleichgewichtes, Manuscript. 

Lesowirz, J. L., and Beramann, P. G., 1957, Irreversible Gibbsian Ensembles, Ann. Phys- 
nee th, ING 

OnsaGcer, L., 1931, Reciprocal Relations in Irreversible Processes, Phys. Rev. 37, 405, 38, 
2265. 

Onsacer, L., and Macuuup, S., 1953, Fluctuations and Irreversible Processes, Phys. Rev. 
91, 1505, 1512. 

SmotucHowski, M., Rirrer von Smoxav, 1908, Molekular-kinetische Theorie der Opalescenz 
im kritischen Zustande, sowie einiger verwandter Erscheinungen, Ann. Phys., 25, 205. 

—— 1916, Drei Vortrage iiber Diffusion, Brownsche Molekularbewegung und Koagulation von 
Kolloidteilchen, Phys. Zeitschr. 17, 557, 585. 

UntHorn, U., 1960 a, Macroscopic observables and generalized canonical ensembles. Arkiv 
for Fysik 17, 233. 

—— 1960 b, On statistical mechanics of non equilibrium phenomena. Arkiv fér Fysik 77, 193. 

Wits, C. R., 1957, Quantum Mechanical Irreversible Gibbsian Ensembles, Thesis, Syracuse. 

Yamamoto, T., 1953, Statistical Mechanics of General Brownian Motions Underlying Irrever- 
sible Processes. Progr. Th. Phys. 10, 11. 

ZERMELO, E., 1900, Anwendung der Wahrscheinlichkeitsrechnung auf dynamische Systeme, 
Phys. Zeitschr. 1, 317. 


Tryckt den 1 april 1960 


Uppsala 1960. Almqvist & Wiksells Boktryckeri AB 


360 


